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Cassiopeia A: dust factory revealed via submillimetre polarimetry 
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ABSTRACT 

If Type-ll supernovae - the evolutionary end points of short-lived, massive stars - produce a 
significant quantity of dust (>0.1 M©) then they can explain the rest-frame far-infrared emis- 
sion seen in galaxies and quasars in the first Gyr of the Universe. Submillimetre (submm) 
observations of the Galactic supernova remnant, Gas A, provided the first observational evi- 
dence for the formation of significant quantities of dust in Type-ll supernovae. In this paper 
we present new data which show that the submm emission from Gas A is polarised at a level 
significantly higher than that of its synchrotron emission. The orientation is consistent with 
that of the magnetic field in Gas A, implying that the polarised submm emission is associated 
with the remnant. No known mechanism would vary the synchrotron polarisation in this way 
and so we attribute the excess polarised submm flux to cold dust within the remnant, pro- 
viding fresh evidence that cosmic dust can form rapidly. This is supported by the presence 
of both polarised and unpolarised dust emission in the north of the remnant, where there is 
no contamination from foreground molecular clouds. The inferred dust polarisation fraction 
is unprecedented (fpoi ^ 30 per cent) which, coupled with the brief timescale available for 
grain alignment (<300yr), suggests that supernova dust differs from that seen in other Galac- 
tic sources (where /poi = 2 — 7 per cent), or that a highly efficient grain alignment process 
must operate in the environment of a supernova remnant. 
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1 INTRODUCTION 

The larg e quantities of dust seen i n high-redshift quasars and 
galaxies jPriddev et all I2003L l2008l : IWang etaLll2007l . l2008h . at 
a time when the Universe was only ~ 1 Gyr old, suggests that a 
rapid mechanism for dust production must operate. Type-Il su- 
pernovae (SNe) are good candidates for these dust factories as 
they evolve to a dust-producing phase in only a few hundred Myr 
and contain a high abundance of heavy elements. Theory pre- 
dicts that each Type-ll SNe shou ld produce ^0.1-1 Mp of dust 
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20071: iBorkowski et al 

20061 : ISugerman et all l2006l : iMeikle et al 



l2006l : IWilliams et al"] 
l2007h . e.g. the maxi- 
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(iTodini & Ferrarall200ll : [Nozawa et al. 2003 ; iBianchi & Schneiderl 
12007) and, if true, this can account for the dust obse rved at high 
redshift ( Morgan & Edmunds 2003 : Dwek etaDl2007h . 

Evidence for this quantity of dust forming in young SNe in 
the local Universe has been scant and controversial. Observations 
from mid-/far-infrared (IR) satellites (IRAS, ISO and Spitzer) detect 
only warm dust and find orders of magnitude less than predicted 
by theory dPwek et al.l Il987l : iLagage et all Il996l : IPouvion et aP 



mum warm dust mass inferred f or Cas A from Spitzer data is 
0.03 - 0.05 M© (iRho et al.ll26o8l) . Submillimetre (submm) obser- 
vations of Cas A seemed to provide the first evidence for large 
(~2M0) quantities of colder dust manufactured in the supernova 
explosion jPunne et al.l |2003[ hereafter D03) . However, subse- 
quent observations of CO and OH towards Cas A suggested that 
some or all of the submm emission may originate from dust in 
a foreground molecular cloud complex r ather than the remnant 
(iKrause et al.ll2004 IWilson & Batrlall2005 ). Although the level of 
foreground contamination is highly uncertain - because of the dif- 
ficulties in converting a molecular line intensity into a submm flux 
density - the molecular data have cast doubt on the idea that signif- 
icant quantities of dust can form rapidly in SNe. 

Another explanation for the submm excess, this tim e with- 
out a significant mass of dusty material, was proposed by IPwekl 
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Elongated iron needles with a much higher submm emis- 
sivity than canonical dust can also account for the far-IR and 
submm spectral energy distribution (SED) of Cas A. If present in 
the general ISM, they should contaminate the polarisation signal 
from the CM B which could have importance consequences (e.g. 
lBowdenetaI[ 2004). 

Polarimetry at 850 /xm is extremely challenging but provides 
one way to test the competing hypotheses for the submm emission 
in Cas A: cold dust in Cas A, foreground contamination, or iron 
needles. We report such observations here, with the paper organ- 
ised as follows: §2 describes the submm and radio polarimetry; §3 
presents the results and investigates the robustness of the submm 
polarimetry; in §4, we use the polarimetry to place constraints on 
the fraction of submm emission arising from the remnant and com- 
ment on the implications for grain alignment theories. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Submm observations 

Our submm observations took place during 2004 October 17-18 
using the SCUBA polarimeter (Murra v et aa]ll997l : [Holland et al.l 
il999 : Greaves et a l. 2003) at the 15-m James Clerk Maxwell Tele- 
scope (JCMT), HawailH We used SCUBA at 850 /xm where the 
array contains 37 bolometers. The instrument comprises a rotating 
quartz half- wave plate and a fixed photo-lithographic-grid analyser 
in a module that is placed in front of the SCUBA cryostat's entrance 
window. The weather conditions were exceptional - very stable 
with very little sky noise, with T23oghz ^ 0.04 throughout. Two 
positions on the bright rim of Cas A were observed (Fig.[T]), one in 
the west (hereafter W) and one in the north (hereafter N). Our orig- 
inal Cas A submm map (D03) was taken in SCUBA's 'scan-map' 
mode; for the polarimetry we were forced to make a 'jiggle' map, 
stepping the secondary mirror of the telescope in a 16-point pattern 
in order to provide a fully sampled image at 850 /xm, whilst chop- 
ping a distance of 180 arcsec at around 7 Hz (the maximum possi- 
ble displacement) - to remove the atmospheric signal, and nodding 
the telescope every 32 s to correct for slowly varying sky gradi- 
ents. The chop positions were chosen to avoid the remnant as much 
as possible; Fig. [T] shows the location of the fields and their chop 
positions. A detailed discussion of the potential effects of chop- 
ping onto emission in the reference positions is presented in §3. A 
full 16-point jiggle was taken at each of 16 positions of the half- 
wave plate, separated by 22.5°. Thus 16 maps complete one full 
rotation, yielding four redundant polarisation observations. We co- 
added five and three such observations for the Western and North- 
ern regions, respectively yielding total integration times of 65 mins 
for the Western region and 45 mins for the Northern region. 

The submm data w ere reduced using the SCUBA User Re- 
duction Facility (SURF -Ijenness & Lightfoot 1998) and POLPACK 
(lBerrv&Gledhillll200lh . The basic analysis consisted of remov- 
ing the nod, flatfielding, correcting for extinction, flagging bad 
bolometers, removing sky noise and clipping noisy data. For the 
extinction correction, we employed a polynomial fit to the Caltech 
Submillimeter Observatory's 230-GHz skydip data. To remove any 



^ The JCMT is operated by the Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom, the 
Netherlands Organisation for Scientific Research, and the National Re- 
search Council of Canada. 




Figure 1. Cas A at 850 /xm, from lD03l showing the regions observed with 
the SCUBA polarimeter. On-source fields are denoted by solid lines while 
chop positions are shown by dashed lines. Position W (white) is the bright- 
est region at radio and submm wavelengths and is coincident with a CO 
emission-line peak. Position N (red) is where the dust-to-synchrotron ratio 
is lowest and where there is no foreground emission. 



residual time-varying sky level, we made preliminary maps in or- 
der to choose bolometers lying in emission-free regions and then 
subtracted their median signal in every second of integration from 
the full array. We experimented with different sky bolometers and 
different thresholds for removing noisy bolometers to assess their 
impact on the polarimetry measurements. We found some small 
changes in the derived polarisation parameters but all maps were 
consistent within the errors (see §3 for more details). 

Instrumental polarisation (IP) arises through the SCUBA 
optics and the GO RE-TEX membrane in front of the dish 
dOreaves et al ]|2003). IP is elevation-dependent and can be char- 
acterised by making observations of unpolarised planets. The IP 
was removed using the SURF task, REMIP, which uses a look-up 
table of bolometer IPs and corrects for the elevation at the time 
of the observation. The uncertainty in the IP removal is ±0.5 per 
cent (Matthews et al. 2001). Once treated in this way, each integra- 
tion was re-gridded separately using a Gaussian weighting function 
with a scale of 7 arcsec (approximately half the FWHM for SCUBA 
at 850 fim) onto 6.2-arcsec pixels (the step size of the jiggle pat- 
tern), ready for the polarisation analysis. From each complete po- 
larimetry observation, four degenerate maps are made for each of 
the four equivalent waveplate angles. The four maps at each an- 
gle are stacked and the variance for each pixel calculated from the 
scatter between the four maps. POLPACK was then used to derive 
the Stokes U, Q and / values for each pixel in each complete po- 
larimetry observation by fitting the following function to the data 
dSparks&Axonll 19991) 

Ik = ^ [/ + e(Ocos20fc + t/sin2(/)fc)] 

where Ik is the expected intensity in image k, (j)k is the position 
angle of the wave plate after correction for the parallactic angle for 
image /c, t is the analyser transmission factor and e is the polarising 
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efficiency factor. The error in each Stokes parameter is derived from 
the input variances. 

These /, Q and U images were then median stacked for repeat 
observations and a final data cube was created. Polarisation frac- 
tions, p, and position angles (PAs), 0, were derived from the Stokes 
parameters (Hildebrand et al. 2000) after correcting for the Ricean 
bias - i.e. the increase in p that results from constraining p to be 
positive even when Q and U are consistent with p = 0. Ricean bias 
is a significant effect for low signal-to-noise measurements (SNR 
<2). 

We define the normalised Stokes parameters and associated 
errors as 



Table 1. Summary of VLA observations. 



1/2 



2 2 2 \ / ^ 



Debiased p values are calculated using 



X 100 per cent 



where the bias is given by 

2 2 I 2 2 



q"^ + u'^ 



(1) 



(2) 



The error in p is ap = a/A x 100 per cent. The polarisation PA - the 
position of the electric field vector measured relative to Celestial 
North - is given by 



^ = 0.5 arctan 



(3) 



The error on is given by (76) = ^ x 28.6° . 

The polarisation data were binned to a resolution of 18 arcsec 
to improve the SNR. Vectors were chosen for the catalogue if 
cTp < 4 per cent, thereby removing unreliable vectors in regions of 
high noise. We have also imposed a cut of SNR ^ 2 in p as vectors 
with below this level have large errors in and t he bias correction 
we m ake is really only applicable at SNR ^ 2 jHildebrand et al.l 
I2OO O). The SNR cut affects only 2/57 vectors and has no conse- 
quences for our results (see ^ 13.21 ). The absolute accuracy in posi- 
tion angle is limited by the systematics of removing the sky noise. 
This was found to be zb6°. 

As discussed bv lGreaves et al.l (l2003h . there is a minimum be- 
lievable polarisation percentage within the main beam of SCUBA 
polarisation maps due to potential contamination of the polarisation 
signal from the extended sidelobes of the JCMT beam. The value 
of the minimum believable polarisation is given by the relation 



Pcrit ^ 2 Psl Psl 



(4) 



where psi is the polarisation in the sidelobe (i.e. off-centre posi- 
tion) and Psl is the power in the sidelobe relative to the main beam 
(measured from planetary observations). The flux at the position of 
an off-centre source relative to the flux at the map centre is given 
by 5'si/5'mb. We have used observations of Uranus, obtained at the 
time of the Cas A observations, to assess the value of pcrit and 
thereby measure the sidelobe polarisation for our science data. 

The highest value of pcrit will be measured for the largest ratio 
of Ssi/Smh- This value arises for our observations of the northern 
rim of Cas A, in which the brightest peak in the field of view does 
not lie at the array centre. For these observations, we measure a 
flux ratio of 1 . 1 relative to the centre of the array. The power in the 



Configuration 


Date 


Bandwidth 


Duration 


A 


2000.9 


6.25 MHz 


12hr 


B 


2001.2, 2001.3 


12.5 MHz 


llhr 


C 


2000.3 


25.0 MHz 


6.0 hr 


D 


2000.7 


25.0 MHz 


4.5 hr 



beam within this annulus (49.0 < Aa'^ + A^^ ^ 61.8 arcsec, 
where a and 6 represent R.A. and Dec.) is 0.014 relative to the main 
beam, and the IP measured within the same annulus is 9.2 ± 4.6 
per cent. The highest value of pcrit is therefore 0.3 di 0.2 per cent. 
We can therefore accept that all values of p in excess of 0.5 per 
cent arise from the remnant itself. Since all measured values exceed 
this threshold by a significant margin, it is clear that none of the 
measurements are due to polarisation in the sidelobe positions. 

2.2 Radio observations 

Cas A was observed with the National Radio Astronomy Obser- 
vatory 's (NRACB) Very Large Array (VLA) in 2000-01 using all 
four configurations from the most extended (A) to the most com- 
pact (D). Data were taken at four frequencies, 4605, 4720, 4860 
and 4995 MHz, as summarised in Table [T] 

Standard calibration procedures were employed, as outlined 
in the AXVS cookbool0 using 3C 48 to set the flux density scale, 
3C 138 to calibrate the polarisation and J2355-I-498 as a local 
phase calibrator. After initial calibration, multiple passes of self- 
calibration were performed to improve the antenna phase solutions. 

The AXVS maximum entropy deconvolution routine VTESS, 
which maximises smoothness in an image, was used to restore the 
total intensity images and the corresponding routine UTESS, which 
maximises emptiness in an image, was used to restore the Stokes 
Q and U images. An 18-arcsec beam (FWHM) was used for decon- 
volution, to match the resolution of the submm data. The standard 
primary beam attenuation correction was applied in the VTESS and 
UTESS routines. The total restored flux density was 788 Jy. Noise- 
corrected, linearly polarised intensity, polarisation percentage, and 
polarisation PA images were made from the Stokes /, Q and U 
maps. The total polarised flux density was 28.9 Jy, resulting in an 
average polarisation percentage of 3.7 per cent. 

To enable a comparison with the submm data, the radio IQU 
datacube was aligned with the submm pixels using the Kernel Ap- 
pUcation Package (KAPPA) task WCSALIGN. Polarisation vectors 
were then calculated in the same way as for the submm data, but 
using bespoke MATLAB scripts. Radio vectors with < 2.4 were 
chosen, providing complete coverage of the regions where we have 
submm vector information. 



3 RESULTS 

The SCUBA 850-/im polarimetry vectors are shown in Fig. [2l 
alongside the 5 -GHz ra dio v ectors. The underlying images are the 
850-/im scan-map from lD03l and the 5-GHz image, both convolved 
to the same resolution as the binned polarimetry data. The length 



^ NRAO is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
3 www.aips.nrao.edu/cook.html 
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850 micron 5-GHz 




Figure 2. Left: 850-/im polarimetry vectors overlaid on the 850-/im map from DOS. Right: 5-GHz polarimetry vectors overlaid on the 5-GHz map. All maps 
are at the same resolution, 18 arcsec. Vectors are rotated by 90° to show the direction of the magnetic field. The orientation of the submm vectors clearly trace 
the approximately radial B field of the remnant. The 5-GHz vectors have been corrected for an average rotation measure of —110 rad . 



of the vectors represents the degree of polarisation, p, and their PA 
shows the direction of the magnetic {B) field (i.e. the E vectors 
measured from the Stokes parameters have been rotated by 90° for 
display purposes). 

These images clearly show that the submm emission is sig- 
nificantly more polarised than the radio emission, with vectors that 
trace the approximately radial magnetic field in the remnant. Both 
traits indicate a source of polarised emission in Cas A above a gen- 
eral extrapolation of the synchrotron polarised flux density. 

We will now investigate in more detail the robustness of this 
result and explore whether the submm polarisation is due to dust or 
synchrotron emission. 

3.1 Robustness of results 

In this section we will describe our investigation into the robustness 
of the data through comparing several reductions of the polarimetry 
data and also through a simulation of the effects of chopping. We 
conclude that the data are robust and that chopping cannot conspire 
to affect our conclusions. Readers who wish to skip the technical 
details can move on to ^ 13.21 

3.1.1 Data-reduction issues 

We explored the effect of using different sky bolometers and clip- 
ping thresholds whilst reducing the submm data The choice of 
sky bolometers is important as several effects can cause a bolome- 
ter's signal to vary with time during a polarimetry observation: 

(i) A changing sky level, which should correlate across all 
bolometers. 

(ii) Bolometers jiggling on and off a bright source. This is the ra- 
tionale behind choosing off- source bolometers to remove sky vari- 
ations across the whole array. 

(iii) Bolometers may lie on regions of the same brightness 
throughout an observation but see polarised emission which 



changes as a function of wave-plate angle, and therefore as a func- 
tion of time within the observation. 

If bolometers displaying effects (ii) or (iii) are used to remove 
sky noise, they could introduce a temporal signal variation to the 
rest of the array and cause systematic changes in the derived polar- 
isation properties. 

We tried two sky-removal techniques: the first used as many 
bolometers as possible, whilst avoiding bright emission regions; 
a second, more conservative approach used only bolometers two 
rows away from bright emission. We found small changes in p 
(Ap = -0.3) and (AO = -7°) but the differences between 
the reductions are consistent with the errors derived from the dis- 
tributions in each map (ap = 1.2, ae = 4°). As the first method 
uses bolometers which may have jiggled onto polarised emission, 
we will use the second, more conservative approach for the rest of 
the analysis. Our overall conclusions are unaffected by this choice 
of sky-removal technique. 

We have a lso compared our maps to the data presented by 
iMatthews et aP I2OO8) where the mean sky level was not removed 
from the array (only the variations are removed) and bolometers 
showing excess noise were flagged aggressively, resulting in the 
removal of around 30 per cent of the data. We have taken a slightly 
different approach, flagging only the four bolometers with clear 
excess noise, plus one that was dead: around 13 per cent of the 
data in total. The tim e series of many bolometers removed by 
iMatthews et al.l (I2OO8 ) displayed momentary spikes but were oth- 
erwise in line with the rest of the array. We de- spiked these by 
flagging any data more than 3a from the mean. Overall, the mean p 
and derived from both reductions are consistent within the errors 
(Ap =1.0 and = —5°). We conclude that the data are robust 
to minor changes in the data reduction procedure. 

3.1.2 Chopping 

We investigated the effects of chopping onto either polarised or un- 
polarised emission, something which can potentially have a serious 
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impact on submm polarimetry data dmidebrand et"aDl200Qh . We 
synthesised Q and U for our reference positions using either the ra- 
dio values scaled to 850 fim or assuming a value for p and for the 
interstellar medium (ISM) surrounding Cas A (these are described 
in more detail below). The Q and U values from our submm po- 
larimetry data were then corrected by adding half of the Q and U 
values found at each reference position. The / values were taken 
from the 850-/im scan-maps, which comprised the iD03i flattened 
and unflattened reductions (hereafter A and B), the former with 
large-scale undulations removed, and the map used by Krauseet 
(2004) (hereafter C). We created a simulated 'chopped' image from 
our measured Q and U values, and a 'chopped' / image from the 
scan map after processing it with the chopping simulator. This is 
intended to represent the data we measure. We then correct for the 
chopping in Q, U and / and produce a 'corrected' map. This as- 
sumes, of course, that the scan-maps provide a good representation 
of unchopped data (they are the best we have). Chopping onto un- 
polarised flux should not affect but will influence p. We have 
thereby simulated the range of plausible structures around Cas A 
to probe the possible systematics. In utilising the unflattened scan- 
maps we have assumed that any structure is real and not an artefact 
of the data reduction. 

In the first instance we assumed that only flux associated with 
the remnant was polarised and used the radio Q and U maps to 
estimate the flux of polarised synchrotron in the reference beams. 
We accounted for a rotation of 24° in going fro m 5-GHz to 8 5 0-/ini , 
in l ine with the rot a tion m easures derived by Ijones et aP (l2003l) 
and [Anderson et a l.^ (1995). We scaled the radio Q and U values 
to 850-/im using the value of the synchrotron spectral index at the 
reference points {a ^ —0.74) and increased the polarised flux in 
the reference positions by a factor three, roughly the average ratio 
of Psubmm/pradio that wc obscrvc, to allow for the reference beams 
being polarised at the same level as the central field. The results 
are summarised in Table [2] We found that chopping onto remnant 
emission polarised at the level seen in the submm would have a 
small effect on the submm vectors, with an average increase in p of 
+0.6 per cent and an average rotation of +0.4°. 

We next investigated the possibility that any foreground 
submm emission in the reference positions is polarised. We allow 
all of the submm emission in the reference positions to be polarised 
at the level found in clouds in the vicinity of Cas A by the Archeops 
submm balloon experiment (Benoit et al. 2004), which is 23 per 
cent at ^ = — 7r/3 and 12 per cent at ^ = 7r/8. It is not certain 
that the Cas A foreground e mission is polarised at this high level 
as the clouds in lBenoit et al. are > 2° away, but we consider this 
as the most conservative case. We find that the difference in the 
chopped and unchopped maps is Ap = 0.3 — 1.2 in field W and 
Ap = 0.2 — 1.9 in field N. The rotation introduced by chopping is 
small, the largest being +1.2° in W and —2.7° in N. These simu- 
lations are shown in Fig. [3] and the results summarised in Tabled 
None of the small changes introduced by chopping alter our con- 
clusions. 

3.2 Submm and radio polarisation properties 

Polarisation vectors were created as described in § 12.11 using the 
same positions in the observed 850-/im and 5-GHz images. We 
also predict the observed polarisation expected at 850 /xm, given 
the measurements in hand at 5 GHz and assuming that the ex- 
cess submm emission is unpolarised. We scale the 5-GHz Q and 
U fluxes using a 1.4 G Hz/5 GHz spectral index map. We use the 
submm scan-map from lD03l as the / image, since this exceeds the 
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Figure 4. Histogram showing the distribution of p for the 5-GHz data 
(black) and the 850-/im data (red). The distributions are significantly dif- 
ferent. The predicted polarisation at 850 /im, using the 5-GHz data as de- 
scribed in the text, is shown by the dashed line. 



predicted synchrotron flux by ~30 per cent (due to dust emission, 
either in the supernova remnant or in front of it). We then calcu- 
late p using the scaled radio Q and U data and the submm / data. 
This is the polarisation we would expect if the submm / emission 
consisted of synchrotron and unpolarised dust. 

A comparison of the p distributions in the radio and submm 
wavebands is shown in Fig.|4]where is it clear that the submm emis- 
sion shows a significantly higher degree of polarisation than does 
the radio. The means of the distributions are Psubmm = 13.5 + 0.8 
per cent and Pradio = 3.7 + 0.2 per cent when selecting vectors at 
p/dp'^ 2 and dp < 4 per cent in the submm waveband. If the SNR 
threshold is ignored, the mean Psubmm drops slightly to 12.9 ± 0.9 
per cent. The means are thus different at the > 8a level. If we com- 
pare to the predicted values at 850 /im (dashed histogram in Fig.|4]) 
the difference becomes even more pronounced. The results of the 
chopping simulations suggest a minor reduction in Psubmm of order 
0.7 per cent may be required, but there appears to be no possibility 
of reconciling the two distributions. 

Having established a clear excess in polarisation at 850 /xm, 
we now wish to understand whether this difference is due to a re- 
duction in the level of synchrotron polarisation between 850 /am 
and 5 GHz via a depolarising mechanism in the remnant, or to 
an extra source of polarised emission at 850 /xm (i.e. aligned dust 
grains). 

Synchrotron emission in supernova remnants should be highly 
polarised (~70 per cent) but observed polarisations at radio fre- 
quencies in young remnants like Cas A are typically much lower, 
implying that the magnetic field is highly disordered on scales be- 
low an arcsec ( Milnelll987i) . Another cause of reduced polarisation 
at radio wavelengths is rotation due to the Faraday effect, in which 
the plane of polarisation is rotated as the polarised waves pass 
through a magnetised plasma. This effect is wavelength-dependent, 
with longer wavelengths suffering more rotation. Emission from re- 
gions at different depths in the plasma of Cas A will be rotated by 
differing amounts as they travel to the observer. This causes 'in- 
ternal depolarisation' at the wavelengths which undergo the most 
rotation. 

Depolarisation at radio frequencies has been studied in Cas A 
and found to be negligible at frequencies of >5 GHz. From 5 GHz 
to 2.2 /im, the fractional synchrotron polarisation remains roughly 
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Al - map A - Archeopsl Q,U 



Al - map B - Archeopsl Q,U 
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Figure 3. Simulated chopped and unchopped images as described in § 13.11 On the left of each sub-panel is the 'chopped' image, meant to represent the 
polarimetry data. Next to it on the right is the corrected image. The upper row is for region W and the lower row is region N. The scan-maps used are A and 
B and the Q and U model is Archeops 1 (see Tableland ^ l3.1| for a detailed description). The effects of chopping on / can be mostly seen at the edges of the 
fields where the edge of the remnant was chopped onto (see Fig.[T). 



Table 2. Results of simulated chopping onto polarised emission. 



W: western field 



N: northern field 



Model 


Submm / map 


Pcor 


Pchop 




AO 


Pcor 


Pchop 


Ap 


AO 


Scaled radio Q, t/ x 3 


A 


11.1 


11.9 


+0.8 


+0.4 


13.0 


13.7 


+0.7 


+0.9 




B 


8.2 


8.7 


+0.4 


+0.4 


11.1 


12.2 


+1.0 


+0.9 




C 


10.6 


11.1 


+0.5 


+0.4 


See note* 








Archeops 1 


A 


10.7 


11.9 


+1.2 


+1.2 


12.5 


13.7 


+1.2 


-0.7 


p = 23 per cent, = — 7r/3 


B 


8.2 


8.7 


+0.5 


+0.9 


10.3 


12.2 


+1.9 


-2.7 




C 


10.4 


11.1 


+0.7 


+0.5 










Archeops 2 


A 


11.6 


11.9 


+0.3 


-0.3 


13.5 


13.7 


+0.2 


+0.5 


p = 12 per cent, ^ = tt/S 


B 


8.5 


8.7 


+0.2 


-0.3 


11.7 


12.2 


+0.5 


+1.6 




C 


10.9 


11.1 


+0.2 


-0.1 










Observed map 


polarimetry 




12.6 ± 1.2 








15.1 + 0.8 







Notes: Column (1) refers to the model used for Q and U in the reference positions; (2) submm scan-map used to determine the reference / value and 
combined with the measured submm Q and U to produce the quoted p values (see text); (3) mean p in the corrected map; (4) mean p in the 'chopped' map, 
error on the mean p is similar to that quoted for the polarimetry in the final row of the table (~ 1 per cent in each field); (5) Ap = Pchop — Pcor; (6) 
AO = 

^chop ~ ^cor- Columns 3-6 are repeated for positions W and N which each have different reference chop positions. * scan-map C cannot be used for 
position N as the background is too negative. 



constant (4-7 per cen t, Ijones et al.ll2003l : lAnderson et al.lll995l : 
iKennev & Dent|[l985b which suggests that the much higher polar- 
isations observed at 850 /xm cannot be due to synchrotron radia- 
tion. In some specialised geometries, depolarisation can lead to an 
increase in fractional polarisation followed by a decrease at suc- 
cessively larger wavelengths, but this occurs only under conditions 
of very strong depolarisation (^ 10 - ICioffi & Jo nes 1980), which 
does not apply in Gas A above 5 GHz. The lack of depolarisation at 



5 GHz thus implies that the high 850-/im polarisations arise from 
an independent source. 

Another piece of evidence for internal depolarisation could 
have come from comparisons of depolarisation with X-ray bright- 
ness, serving as a pr oxy for the density of the depolarising ther- 
mal electrons. Indeed. lAnderson et al.l (4995) showed a very strong 
correlation between the X-ray brightness and the depolarisation of 
emission between 5 GHz and 1.4 GHz, arguing for internal depolar- 
isation at the low radio frequencies. However, no such dependence 
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is seen in our comparison of depolarisation between 850 /xm and 
5 GHz and the X-ray brightness (Fig. [5] left). Thus, the high frac- 
tional polarisations at 850 \xm must arise from another cause. This 
850-/im-independence is also demonstrated in the plot of Pradio vs. 
Psubmm in Fig. [5] (right) where no correlation is seen. 

We now turn to the PA of the polarisation. In Fig. [Gfleft) we 
show the PAs (0) of the radio versus the submm. Filled symbols 
indicate the Western region and open symbols the Northern re- 
gion. There is a very good correlation between PAs contrary to 
the lack of correlation between polarisation strengths. This sug- 
gests that the submm polarisation is largely arising from a different 
physical mechanism than the synchrotron radio polarisation. How- 
ever, both polarisation signals reflect the same underlying magnetic 
field in Cas A. This is further demonstrated in Fig. [6[right) where 
we shows the PAs of the rad io and submm and 2.2- fim polarisa- 
tion vectors ( Jones et al as a function of azimuth. The radio 
measurements have been cor rected for the local Far aday rotation, 
redetermined from the data of Anderson et al. ('1995'). They quoted 
an average value for the rotation measure (RM) of —110 rad m~^, 
while we determine a range of values from —50 to —302 rad 
for our range of azimuths. The extrapolated rotations at 850 /xm are 
<0.01 degrees, so no correction for Faraday rotation has been made 
for those, or for the 2.2-/im data. 

There is a very good correlation between all three sets of mea- 
surements, suggesting that whatever is responsible for the extra 
polarisation signal in the submm waveband is related to the same 
magnetic fields that are responsible for the synchrotron emission at 
5 GHz and 2.2 /xm and, perforce, is intrinsic to Cas A and is not 
due to foreground material. The small mean offset between the cor- 
rected 5-GHz PAs and those at 850 /xm and 2.2 /xm are consistent 
with the systematic uncertainty in the RM correction. To first order 
the PAs all provide strong evidence for a radial magnetic field , as is 
common in young supernova remnants ( Anderson et al.l l 19951) . 

If we take the submm and 2.2-/im magnetic field PAs at their 
face value, ignoring the systematic uncertainties, they agree very 
well with each other and are systematically offset from a radial 
angle by ^10-20^ in a clockwise direction. At present, there is no 
apparent reason for this offset. 



4 DISCUSSION 

Since we cannot explain the submm polarisation as arising from 
synchrotron emission we will now consider the most likely alter- 
native - that it is due to dust aligned with the magnetic field in the 
remnant. We will estimate pa, the polarisation fraction of the dust 
in the remnant. We know that the submm flux we measure is a com- 
bination of radio synchrotron and thermal emission from dust, such 
that /s = /r + /d- Similarly we can write that Qs = Qr + Qd and 
Us = Ur -\- Ud- We can also express Q and U in terms of /, p 
and ^, as Q = Ip cos{20) and U = /psin(2^). Combining these 
equations and solving for dust polarisation, pd, gives: 

[(/sPs - IrPrf + 2/sPs/rPr(l - COs2(6's - Or))Y^^ 

Pd (5) 



r 



n 



20 



40 

polarization % 



60 



80 



which for = Or simplifies to 



Figure 7. Histogram of dust polarisation fraction, pd, calculated using 
Eqn.|6] 



of Pd is shown in Fig. |7] Eqn. [6] assumes that the polarised syn- 
chrotron and dust radiation is emitted with the same intrinsic PA 
(i.e. aligns with the B field in the same way) and that any rotation 
between 5-GHz and 850-/im vectors is due to Faraday rotation at 
5 GHz. This is a reasonable assumption given Fig. [6] We obtain a 
slightly higher value for pd (33 per cent) if we use the difference 
between the RM corrected 5-GHz PA and the submm PA in Eqn.|5] 
instead. 

We can create a map of 850-/im dust flux density, /d, by scal- 
ing the 5-GHz map to 850 /xm as described in ^ 13.21 Subtracting the 
scaled synchrotron map from the 850-/im image leaves the 850-/im 
emission due to dust. This dust maf|j is shown in Fig. [S] with the 
dust polarisati on calculate d from Eqn. (6] Over-plotted are C0(2-l) 
contours from lEales et al.l (in preparation). This shows three things: 
first, there is dust emission in the northern part of the remnant, but 
no corresponding CO emission; second, the level of submm polar- 
isation anti-correlates with the CO contours; the dust polarisation 
be ing lowest at the w ester n peak, consistent with the suggestion 
bv lKrause et"aD (l2004l) and lWilson & Batrlal (l2005h that its submm 
emission is contaminated by foreground clouds; third, there is po- 
larised dust emission in both the north and the west. 

The polarisation signal is correlated with the magnetic field in 
the remnant, which means that some of the emission in the western 
region must be associated with the remnant. None of dust intensity 
or polarised flux in the north is associated with molecular material. 
Thus the finding of Kraus e et al.l (|2004|) that there is 'no cold dust 
in Cas A is demonstrably incorrect. The average dust polarisation 
around the western peak is 13.7 ± 1.4 per cent; in the north, away 
from the CO, it is 39.4 zb 3.2 per cent. If we assume that the in- 
trinsic value of Pd is similar throughout the remnant, we can use 
the uncontaminated values in the north to estimate the contribution 
of the foreground material to the western peak. Since p oc 1 //, 
we simply correct the dust flux in the western peak by the ratio 
13.7/39.4 = 0.34, which gives a flux density of ^ 0.8 Jy intrin- 
sic to Cas A at the Western peak. In order to make a comprehen- 
sive correction for the intervening foreground material a full set of 
molecular gas tracers must be observed over the whole remnant. 



Pd = {Ps -P^)j^ 



(6) 



The fractional dust polarisation calculated using Eqn.[6]is un- 
precedented, with an average pd = 30 it 2 per cent. A histogram 



This map differs slightly from that in as we have used a more ac- 
curate, spatially varying spectral index to extrapolate the radio synchrotron 
flux to the submm waveband. 
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Figure 5. Left: Depolarisation from 5 GHz to 850 fim as a f unction of X-r ay brightness. There is no clear relationship, unlike the 1.4-GHz to 5-GHz depolarisa- 
tion which correlates strongly with X-ray brightness ( Anderson et al. 1995). Right: Radio polarisation versus submm polarisation. The line shows a one-to-one 
correspondence. Filled black diamonds are the Western region while red open circles are the Northen region. There is no obvious correlation. 
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Figure 6. Left: Radio versus submm position angles. Filled symbols represent the Western region while open symbols represent the Northern region. There is a 
very good correlation between the direction of the radio and submm polarisation vectors. Rig ht: Inferred magnetic fie ld angle (PA— 90°) versus azimuth angle 
(0° for celestial north). Solid diamonds are 850-jun i PAs; red triangle s are 5-GHz PAs from [Anderson et al.l Il995h corrected for the local RM as described 
in the text; blue open circles are 2.2-/im PAs from I Jones et al.l (l2003h . Systematic errors bars are shown in the top left. The solid line shows a pure radial 
alignment. 



This data set does not yet exist and therefore we cannot make a 
definitive measurement of the total submm flux associated with 
the remnant at this time. We can, however, estimate a conserva- 
tive lower limit by assuming that only the polarised flux density at 
850 iim{r^ 30 percent) is intrinsic to the remnant. The integrated 
flux density from the dust map in Fig. [8] is 5*850 = 20.1 m and 
taking the average 30 per cent polarised fraction as the minimum 
gives us a total lower limit on the dust flux for Cas A at 850 /xm of 
6.0 Jy, based on our submm polarimetry. 

To convert this flux density into a dust mass we must assume 
a value for the dust mass absorption coefficient a^sso • In our pre- 
vious work (D03) we used a value which was derived from labo- 
ratory studies of cosmic dust analogues which were amorphous, 
non-spherical or aggregates. This gave a high value of a^sso = 

^ This is higher than the value of 15.9 Jy quoted in lD03l due to the lower 
synchrotron contribution at 850 /im when using a varying radio spectral 
index 



0.76 kg If we apply this value to our flux measurement and 
use a temperature of = 20 K (which comes from fitting the 
SED from the IR-submm) we get a dust mass estimate for Cas A of 
~ 1.0 M0. This value is consistent with theoretical prediction s for 
supernova dust yield s (e.g. Todini & Ferrara 2001; Nozawa et alj 
l2003l : lBianchi & Schn eider 200 71) and sufficient to explain th e dust 
at high redshift d Morgan & Edmundsl2003l : lDwek et al.l2007h . This 
is also an upper limit to the possible mass of dust in Cas A given 
the masses of condensible elements formed in supernovae with pro- 
genitor masses in the range applicable to Cas A (13-20 M©- as 
sugge sted by the recent d etermination of the type of Cas A as a 
Ilb bvlKrau se etal.' ('2008')). According to nucleosynthetic models 
bv lLimong i & Chieffi ( 2003), supernovae with progenitors in this 
mass range could produce 0.5-1.0 MqoI dust if the condensation 
efficiency were 100 per cent. 

We can also look at new values for a^sso which have 
been derived from models of dust formation in supernovae sim- 
ilar to Cas A. Values range from /^sso = 0.2m^kg~^ from 
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35 23:23:30 25 
RA 

Figure 8. Synchrotron- subtracted 850-/im image, showing the location 
of the dust. Overlaid are the dust polarisation vectors calculated using 
Eqn.[6] Also show n in blue are CO(2-l) contours, at the same resolution 
teales et al .l2008L (in preparation)). CO contours coincide with the western 
dust peak, but the dust polarisation at these positions indicates that some of 
the 850-/im emission must be due to dust in the remnant. The decrease in 
dust polarisation at the CO peak in the west is consistent with some fraction 
of the submm emission at this location being due to foreground material, 
associated with the CO. Note that there is no CO emission coincident with 
the dust emission in the north. 



iBianchi & Schneided fcOOVh to /^ssq = 0.049 m^kg"^ fo r a 
Cas A-like case from ('Nozawa e t al.ll2003l : iKozasa e t"aDl2008l . in 
preparation). These n values produce dust masses in the range 3.8- 
15 M0. Given the above discussion on the available mass of con- 
densible elements, these values are not physically plausible. This 
suggests that the dust which is emitting the polarised radiation at 
submm wavelengths is not currently predicted by the dust forma- 
tion models and is not necessarily the same as the dust which is 
emitting at the shorter FIR wavelengths seen by Spitzer and IRAS. 
The supernovae formation models currently only deal with spher- 
ical grains, but in order to produce polarised emission the grains 
must be elongated to some degree. Further development of the dust 
formation models may be required in order to produce grains which 
can reproduce the quantity and polarised nature of the submm emis- 
sion in Cas A. In order to have a physically plausible mass of dust 
the emissivity of the dust must be high - at least as high as that 
found in the laboratory studies of amorphous and aggregate grains. 

The highly elongated iron needle grains proposed by IPwekl 
alleviate the emissivity issue, as for long enough axial ratios they 
have emissivities many orders of magnitude higher than spheroidal 
non-conducting grains. Graphite whiskers may also provide a sim- 
ilar solution, especially in light of their recent discovery in mete- 
orites (Fries & Steelell2008b . Iron or graphite needles could reduce 
the dust mass required although a full analysis of the properties of 
such grains which match the new observations is beyond the scope 
of this paper. We therefore do not rule out such exot ic grains and 
will c onsider them in more detail in a future paper (iGomez et al.l 
I2OO8I . in preparation). 



4.1 Grain alignment 

What does this polarisation fraction tell us about the proper- 
ties of the dust in Cas A? Average polarisation fractions in 
typical inter stellar and molecular c l ouds are of the order 2- 
7 per cent (iHildebrand et alJ |200QI : iMatthews & WilsonI l2002l : 
ICurran & Chrvsostomou 2007) though some clouds have values as 
high as 15-20 per cent ( Benoit et al. 2004). We are unaware of any 
measurements as high as those quoted here, which suggests that ei- 
ther the dust or the alignment mechanisms in Cas A differ - perhaps 
unsurprisingly - from those in the general ISM. 

Observations of polarised starlight and thermal emission from 
dust mean that there must be a population of dust grains which are 
non- spheroidal and a mechanism capable of aligning the grains so 
that an appreciable polarisation signal can be detected. The details 
of how interstellar dust aligns so efficiently, despite impacts by gas 
atoms acting to randomise their orientations, h as kept theorists busy 
for half a ce ntury. The original proposal by IPavis & GreensteinI 
(I1950IIT951I) suggested paramagnetic alignment, where the axis of 
rotation of the slightly elongated grains becomes aligned with the 
B field though magnetic dissipation. This theory was unable to ex- 
plain the polarisation levels observed in the ISM and star-forming 
regions because the timescales for alignment are longer than the 
timescales for randomisation by collisions with the gas. 

There have been many significant improve ments to the the- 
ory of grain alignment in recent years - see iLazarianl (l2007h 
for an excellent review. Several mechanisms are now thought 
to operate in different astrophysical environments. Briefly these 
comp r ise radiative torque s ( Dolginov 1972; Draine & Weinsartner] 
Il996l : iLazarian & Hoand Eo07al : iHoang & LazarianI l2008al), me- 
chanical alignment in a s upersonic or subsonic gas flow (iGoldl 
1952; Roberse et al. 1995; ILazarianl [19971 : lYan & LazarianI l2003l : 
Lazarian & Hoang 2007b) and super-paramagnetic alignment 
(Lazarian & Hoang 2008). 

The strong and turbule nt magnetic field in Cas A {B ^ 
0.5mG -'Wrigh t etall ll999^. together with the large abundances 
of heavy elements and the presence of a hot X-ray plasma, means 
that certain mechanisms may be particularly applicable to the en- 
vironment of Cas A. Pinwheel torques arising from the interactio n 
of grains and electrons in a hot plasma (iHoang & Lazarianll2008bh . 
supersonic mechanical alignment, sub- sonic mechanical alignment 
through MHD turbulence ( Yan & Lazarian 2003) or grains with 
super-paramagnetic (e.g. Fe) inclusions jLazarian & Hoand l2008h 
can all lead to highly efficient alignment. The timescales they have 
to operate is small in astrophysical terms, as the explosion of Cas A 
occurred only ~300 years ago. We note that the low synchrotron 
polarised fraction indicates a turbulent and disordered magnetic 
field on sub-arcsecond angular scales. This acts to reduce the ob- 
served polarisation as random orientations are averaged within the 
synthesised beam. If the dust grains were aligned with this tangled 
small-scale magnetic field then we might expect a similar beam de- 
polarisation effect to operate at submm wavelengths. However, this 
would give rise to unphysical intrinsic dust polarisation fractions 
and so we believe that the alignment mechanism must be operating 
on larger scales. The radial B field in Cas A is ordered on large 
scales and thought to arise fr om Ra leigh-Taylor instabilities at the 
reverse shock boundary (Gullll 1975b . These large-scale radial mo- 
tions may also be responsible for aligning the dust in Cas A. 

Even if full alignment is theoretically possible, the line of sight 
angle to the B -field in Cas A together with the random alignment 
of orientations along the line of sight and within the beam should 
all act to reduce the measured polarisation. That we still observe 
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such a high fpoi indicates that the grains themselves are emitting 
strongly polarised radiation. For ideal alignment and line-of- sight 
viewing conditions fpoi ^ 30 percent can be acheived with only 
moderately non- spherical grains of astronomical silicate, with axial 
ratios of between 1.4-1.7 for oblate and prolate grains respectively 
jPadoan e t al. 2001). Thus the level of polarisation per se is not, on 
its own, strong evidence for needle like grains. An estimate of the 
aligned fraction would be required in order to constrain the shape 
of the grains. 

A full application of the various grain alignment models to 
Cas A is beyond the scope of this paper but these data should pro- 
vide an interesting test of aligment theories and in time should also 
shed light on the nature of the grains responsible for the polarised 
submm emission. 



SUMMARY 

• We have discovered unprecedented levels of submm polarisa- 
tion towards the Cas A supernova remnant, significantly in excess 
of that expected from the radio synchrotron, and correlated with the 
magnetic field direction in Cas A. 

• There is no currently known way to produce the polarised 
submm emission from a synchrotron process; a depolarising mech- 
anism capable of producing a peak in /poi as a function of wave- 
length would be required. 

• The polarised emission is therefore strong evidence that a sig- 
nificant fraction of the cold dust detected by DOS is associated 
with the supernova remnant. Further strong evidence is that there 
is both polarised and unpolarised dust emission in the northern re- 
gion which is not affected by foreground molecular clouds. 

• A high dust emissivity is required in order to stay within the 
constraints of the mass of condensible elements available to form 
dust. This limit on the emissivity is not consistent with that cur- 
rently predicted for dust formed in supernovae. 

• The strength of the polarisation signal is unprecedented in 
the general ISM and indicates a highly efficient alignment mech- 
anism at work in Cas A. These data provide a stringent test of grain 
alignment mechanisms, given the very short timescale available, 
-300 yr. 

• Alternatively, the strong signal may arise in exotic grains such 
as iron or graphite needles. Such grains would naturally produce a 
higher emissivity and keep within the mass budget. 

• Higher-resolution polarimetry at 850 /xm is required to con- 
firm the polarisation signal and provide a comparison at the spa- 
tial resolution of the radio synchrotron emission. Measurements 
of the high-frequency (100-200-GHz) synchrotron polarisation are 
required to test whether a perverse radio depolarisation mechanism 
is at work. 

• Future far- IR and submm polar imeters, suc h as those planned 
for SCUBA-2 (Bastien et al.l 12005 ') and FIRI (Hel mich & IvisonI 
^08), will be able to image the submm polarisation in Cas A at 
higher resolutions and frequencies. 
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